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Experimental studies of magnetoresistance in thin superconducting strips subject to a perpendic- 
ular magnetic field B exhibit a multitude of transitions, from superconductor to insulator and vice 
versa alternately. Motivated by this observation, we study a theoretical model for the transport 
properties of a ladder-like superconducting device close to a superconductor-insulator transition. 
In this regime, strong quantum fluctuations dominate the dynamics of the vortex chain forming 
along the device. Utilizing a mapping of the vortex system at low energies to one-dimensional (ID) 
Fermions at a chemical potential dictated by B, we find that a quantum phase transition of the 
Ising type occurs at critical values of the vortex filling, from a superconducting phase near integer 
filling to an insulator near 1/2 -filling. The current-voltage (7 — V) characteristics of the weakly 
disordered device in the presence of a d.c. current bias I is evaluated, and investigated as a function 
of B, I, the temperature T and the disorder strength. In the Ohmic regime (J/e <C T), the resulting 
magnetoresistance R(B) exhibits oscillations similar to the experimental observation. More gener- 
ally, we find that the I — V characteristics of the system manifests a dramatically distinct behavior 
in the superconducting and insulating regimes. 



I. INTRODUCTION 

In superconducting (SC) systems of reduced dimen- 
sionality (i.e., thin films and wires), transport properties 
are strongly affected by fluctuations in the superconduct- 
ing order parameter. The most prominent manifestation 
of the role of fluctuations is the appearance of a finite dis- 
sipative resistance below the mean-field critical temper- 
ature T c of the bulk superconductor. This failure of the 
hallmark of superconductivity - the zero-resistance char- 
acter - may persist to very low temperatures T <C T c , 
where pair breaking is negligible and the electronic state 
can still be described in terms of complex order parame- 
ter field representing the Bosonic degrees of freedom. In 
this regime, while fluctuations in the amplitude of the or- 
der parameter are suppressed, fluctuations in the phase 
field play a dominant role. In particular, when topolog- 
ical defects (vortices and phase-slips) develop dynamics, 
a dissipative voltage is generated in response to a current 
bias. In the T — > limit, their quantum dynamics dom- 
inates and may lead to the formation of a liquid phase, 
characterized by a m etallic or insulating behavior of the 
electronic systerrfS^]. 

In the one-dimensional (ID) case, i.e. in SC wires of 
width and thickness smaller than the coherence length £, 
the resistance essentially never vani she s at finite T due 
to thermal activation of phase-slip s^^ ( for T < T c ) or 
their quantum tunnelling at lower T^isM j n contrast, in 
the two-dimensional (2D) case (SC films), superconduc- 
tivity is well-established at sufficiently low T. However, 
by tuning an external parameter which leads to prolifer- 
ation of free vortices, it is possible to drive a quan tum 
(T — > 0) superconductor-insulator transition (SIT)MI. 
Employing the concept of charge-flux dualitj^, one may 
relate the conduction properties of the electronic system 
to the various phases of vortex matter by interchanging 
the roles of current and voltage. Thus the SC phase is 



associated with a vortex solid, while the insulator can be 
viewed as a vortex superfluid. 

Experimentally, one of the most convenient ways to 
induce a tunable SIT in SC films is by application of a 
perpendicular magnetic field B. At fixed T, a positive 
magnetoresistance R{B) is typically observed in a wide 
range of B. The SIT is then identified in the data as a 
crossing point of these isotherms at a critical field B c , 
separating a SC phase for B < B c from an insulating 
phase for B > B c . At finite T, in both phases the resis- 
tance is typically finite, and the distinction between the 
phases is deduced from the trend of R vs. T: dR/dT > 
indicates a superconductor, and dR/dT < an insulating 
behavior. 

Recent experimental studies of InO devices character- 
ized by a strip geometry^ - namely, a SC wire of width 
comparable to £ - offer an opportunity to probe the 
crossover from a ID to 2D quantum dynamics of the topo- 
logical phase-defects. The prominent observation is that 
in the presence of a perpendicular field B, the magne- 
toresistance R(B) exhibits oscillations which amplitude 
is sharply increasing at low T, in striking resemblance 
to the behavior of Josephson arrays^ and SC network 
systems^. Moreover, the SIT at a high field B c appears 
to be preempted by a multitude of transitions at lower 
fields, from a SC to an insulator or vice versa alternately. 
These are indicated by multiple crossing points between 
different isotherms R(B). 

The periodicity of the above mentioned oscillations is 
consistent with a single flux penetration to the sample. 
This suggests that the observed SC or insulating behavior 
of the system is determined by commensuration of vor- 
tices within the strip area. In particular, when an integer 
number of vortices can be fitted along the strip form- 
ing a uniformly-spaced chain, superconductivity may be 
supported even at sufficiently high B such that a large 
fraction of the sample area turns normal. However, devi- 
ation from commensurability of the vortex filling forces 
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FIG. 1: (color online) (a) Top view on a chain of vortices 
in a superconducting strip, (b) The line-junction model for 
the system; purple dashed lines represent the Josephson cou- 
pling gj, and brown dotted lines the Coulomb interaction U 
between the two SC wires. 



a frustrated vortex configuration, thereby weakening su- 
perconductivity. In this case, the quantum mechanical 
character of vortices is manifested by the formation of de- 
localized vortex states, facilitating their mobility across 
the width of the stripP*] As a consequence, the tuning of 
vortex filling away from commensurability can possibly 
induce a quantum phase transition to a liquid state, of a 
metallic^ or insulating character. This commensurate- 
incommensurate effect may also be manifested as magne- 
tization plateaux, as was predicted in a theoretical study 
of bosonic ladder^. 

In a recent papei^ we have studied this phenomenon 
within a theoretical model for a quantum vortex chain in 
a ladder-like SC device (see Fig. [T]), which particularly 
addresses the strongly quantum fluctuation regime where 
the parameters are close to a SIT. It was shown that such 
system may exhibit multiple quantum phase transitions 
of the Ising type, manifested as SC-insulator oscillations 
of the Ohmic resistance R(T,B). This reflects an inti- 
mate correspondence between charge-flux duality across 
a SIT, and the order-disorder duality characterizing the 
Ising transition at 1 + 1-dimensions. 

In this paper we present a detailed theory for the elec- 
tric transport properties of the quantum vortex chain in a 
weakly disordered SC ladder. In particular, we derive the 
current -voltage (I — V) characteristics of the device in 
the presence of a d.c. current bias I, and investigate their 
behavior as a function of B, /, the temperature T and the 
disorder strength. We find that the I — V characteristics 
of the system manifests a dramatically distinct behavior 
in the SC and insulating regimes. In the Ohmic regime 
(7/e <C T), this yields an oscillatory magnetoresistancc 
i?(T, B) which exhibits T-dependence compatible with 
the experimental data. 

The paper is organized as follows: in Sec. In] we con- 



struct the line-junction model for the SC strip, and de- 
rive its mapping to ID Fermions and consequently to 
the quantum Ising chain. In Sec. |III| we provide a de- 
tailed calculation of the dissipative voltage in a current- 
biased strip, and derive expressions for the non-linear 
I— V characteristics and T-dependent magnetoresistance 
in the various regimes (the SC phases, insulating phases 
and critical regions). Our conclusions and discussion of 
the relation to further experiments are summarized in 
Sec. El 



II. THE MODEL 

We consider a thin SC strip of length L>( and width 
w > £, subject to a strong perpendicular magnetic field 
below the 2D SIT (i.e., B < B c ). A ID chain of vortices 
is formed along the central axis of the strip, which can 
be viewed as a ID system of particles in the presence of 
a self-organized effective potential dictated by the com- 
bination of vortex-vortex interaction and the boundary 
conditions [Fig. [lja)]. In particular, the interface with 
the vacuum at the strip edges induces an effective " image 
charges" potential, and bulk-superconductor contacts 
connected to both ends of the strip enforce a fixed phase 
of the SC order parameter at x = ±L/2. As a result, 
the effective potential acquires the form of a periodic ID 
lattice of pinning sites separated by a uniform spacing 
a = L/N, where N = l[BwL/$ ] (with $ = hc/2e the 
flux quantum, and I[z] the integer value of z) denotes 
the total number of vortices^!. Assuming further that 
the high vortex density in this case leads to near merg- 
ing of their cores along the central axis of the strip, the 
system becomes essentially equivalent to a line-junction 
formed by a pair of parallel SC wires separated by a nor- 
mal barrier [Fig. fTKb)], subject to a magnetic field B 
perpendicular to the junction plane. 

In the low T regime, pair-breaking is negligible and 
the properties of this system are dominated by quantum 
phase-fluctuations of the SC condensate. It is therefore 
possible to model it as a 2-leg bosonic laddei^ (or, equiv- 
alently, a ladder-like Josephson array^^) , where a coordi- 
nate x = ja (j integer) denotes the locations of vortex 
cores in the continuum limit. The dynamics of the col- 
lective phase field in the wires (<j> n (x, t) with n = 1, 2) is 
governed by the effective ID Hamiltonian 



Ho — H i + H 2 + Hi n t , 
in which (using units where h = 1) 



H n 



i n 



dx 



U p 2 n+^(d^n) 2 



(1) 



(2) 



H^t = J ^dx {-gj cos(0i - 4> 2 - qx) + U pip 2 ] .(3) 
Here the operator p n {x) denotes density fluctuations of 
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Cooper pairs in wire n, and can be represented as^ 

Pn(x) = --d x e n {x) +p Ve' 2 **- 9 "' (4) 

7T * — ' 

in terms of the conjugate field 9 n (x) satisfying 
[4> n (x), d x 9 n (x')} = i , k5{x' — x). The first term in Eq. 
([2]) hence describes a charging energy; p s is the super- 
fluid density (per unit length) assumed to be monotoni- 
cally suppressed by increasing B, po — p s (B — 0) and m 
is the electron mass. The inter-wire coupling [Eq. 
consists of a Josephson term and an inter-wire Coulomb 
interaction, of coupling strengths gj and U, respectively. 
Finally, the parameter 
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(5) 



parametrizes the deviation of the vortex density from the 
closest commensurate value, i.e., it denotes vortex "dop- 
ing". We note that H describes an ideal system, to 
which we later add a disorder potential. 

To further analyze the properties of this model, it is 
convenient to introduce symmetric and antisymmetric 
phase and charge fields via the canonical transformation 



^ = 4?(0i±02), e ± = ± w ■ ( 6 ) 



In terms of these variables, the Hamiltonian ([!]) is sepa- 
rable: 



H = H + + H_ 



where 



H + =H 



(+) 

LL > 



H- = H 



(-) 

LL 



dx 



-gj cosCv^- - qx) + g c 



cos 



) ; (7) 



H 



(±) 

LL 



V± 

2tt 



dx 



K ± {d x e ± ) + —{d x <t> ± ) 2 



(8) 



This yields 



and the parameters are given by 



K+ = 



Um(U ± U) 



K 2 Ps 



V± 



Ps(U Q ±U) 



Am 



9c 



2Up 2 



(9) 
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[{d x ef + (d x( f>) 2 



(ii) 



+ I dx [—gj cos(2(/) — qx) + g c cos(26*)] 



Here we have accounted for the most relevant interaction 
terms, neglecting umklapp terms included in the charging 
energy which are effectively suppressed due to the rapidly 
oscillating factor in Eq. Q. The symmetric mode (cor- 
responding to the plasmons of total charge) governed by 
H + is therefore gapless. However, the behavior of the 
antisymmetric mode is dictated by the competition be- 
tween two interacting (cosine) terms, and depends cru- 
cially on the value of the Luttinger parameter . Below 
we focus on the regime of parameters close to a SIT in 
ID wires, where quantum fluctuations in the phase and 
charge fields are maximized; i.e., A'_ ss K c = 2 (see Ref. 
©. 

Wc next define new canonical fields 



1 

r(") 



6 = 



(10) 



in terms of which Hj^ L acquires the form of a Luttinger 
Hamiltonian with an effective Luttinger parameter K — 
K-/2. For K- close to K c = 2, we thus obtain K « 1. 



This model can be refermionized by introducing right (R) 
and left (L) moving spinless Fermion fields^! 



±ik F x p i(T<i>+6) 



(12) 



in terms of which _H_ becomes a free Hamiltonian. Here 
the short-distance cutoff a is set by the lattice constant a 
characterizing the vortex chain, and the "Fermi momen- 
tum" kp — n/a + q is determined by the vortex filling 
factor [see Eq. dfj])]. Quite interestingly, this implies that 
near a SIT, it is natural to adapt a duel representation 
of this system in terms of fermionic vortex fields. This 
stems from the approximate self-duality of (i.e., its 
symmetry to exchange of 4> and 9), implying that the nat- 
ural degrees of freedom are composites of a pair charge 
(2e) and a unit of flux quantum. 
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The fermionic representation of -£/_ is given by 

H- = J dx{v-{ip R (x)(~id x )ip R {x)-ipl(x)(-id x )ip L {x)} 

- Vv[iPr(x)iPr(x) + ipl(x)i/j L (x)} 

- J[^ R {x)i) L {x) + ^ L (x)il) R {x)\ 

+ V[^ R {x)^ L {x)+^ L {x)^ R {x)]} (13) 

where J — nagj, V — irag c and the vortex chemical 
potential is [i v — nv-q 7 which vanishes at commensu- 
rate fillings. Following the analogous problem of spin- 
1/2 ladders^E^, it is useful to decompose the complex 
Fermions [Eq. ( 12 )] in terms of the Majorana fields 



(lp v + 1(>1) , f]2v = ~~J^ ty» - 4) ( 14 ) 



[y = R,L). Recasting Eq. (13) in fc-space and using the 



Fourier transformed fields T)j v ,k — f}j v we obtain 



I v_k iA ( ° ] 



-iA 



(o) 



ifj, v 

V-k — ifi v 
-ifi v u_fc —iA d 
V i[j, v iA { d 0) / 

( WlR,k , V2L,k , r\2R,k , ViL,k ) 



here 



(15) 
(16) 



denote the gaps in the excitation spectrum for commen- 
surate vortex filling (/j, v — 0), in which case "Hfc decouples 
into two independent blocks. Since J, V are positive, the 
u sector is higher in energy. 

We now focus on the case of interest, where the system 
is assumed to be in the SC phase but close to a SIT so 
that the Josephson energy J is slightly larger than V, and 
A^ <C Au \ In this case, the high energy sector u can 
be truncated, and the low-energy properties are governed 
by the cJ-type Fermions. Most notably, the gap A d °^ can 
change sign upon tuning of J below the critical value 



J r = V where A 



(o) _ 



0. Indeed, for fi v — each species 



of free massive Fermion models described by ( 15 ) can be 



independently mapped to an Ising chain in a transverse 
fielcPED. j n particular, the low energy sector d can be 
described by the spin Hamiltonian 



2^ J 



(17) 



which possesses a quantum critical point at J = V. 

When finite vortex "doping" is introduced by tuning 
B away from B^ such that fx v ^ Q, the original d and 
u sectors mix. However, the resulting long wave-length 



theory can still be cast in terms of two decoupled sec- 
tors denoted d (low) and u (high). Moreover, the energy 
spectrum 



e„,d(fc) 



} 2 _k 2 ±2 



J 2 V 2 + (fi v v-) 2 k 2 



reduces in the k — > limit to the same form as the fi v 
case: 



1/2 

(18) 

= 



£ u ,d{k) ~ A Utd + 



1 <d fc2 

2 A u h 



with the modified velocities 



u u,d 



1± 



fly 
JV 



and modified gaps given by 

A u , d (B) = J±V . 



(19) 



(20) 



(21) 



The B-dependence of A Uj d is oscillatory due to the de- 
pendence of V on the vortex doping fi v [Eq. ( 18 )]. While 



A u remains positive and large for arbitrary fj, v , a quan- 
tum phase transition occurs at a critical value of fi v 
[which can be traced back to a sequence of critical fields 

B^ via fJ, v (q) and Eq. (|5j)], where changes sign. As 



B -> B, 



(N) 



, one expects the scaling 
\AA ~ \B~bW\ 



(22) 



As we show in the next Section, the above discussed Ising 
like quantum critical points correspond to SC-insulator 
transitions, marked by a dramatic change in the trans- 
port properties. 



III. I-V CHARACTARISTICS AND 
MAGNETORESISTANCE 

We next study the transport properties of the system 
in the presence of a weak scattering potential, generi- 
cally induced by random, uncorrelated impurities along 
the coupled wires. To this end, we include a linear cou- 
pling of the density operator p n (x) [Eq. Q] to a disorder 
potential Vd(x) in the Hamiltonian. The leading contri- 
bution to dissipation arises from the backscattering term 
of the forirP 



H D = 



71=1,2 



dxQ n (x) cos{29 n {x)} 



(23) 



where we assume 

(C»(a:)) = , (Cn(z)Cn'(z')) = D5(x - x')6 n , n ,. (24) 

Here and throughout the rest of the section, the definition 
of ( ) includes disorder averaging. As a result of phase- 
slips generated by Hjj , a finite voltage will develop along 
the SC strip when driven by a current bias /. 
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To introduce a d.c. current bias I, we add a time- 
dependent term It to the total charge operator 

Q = --(0i + o 2 ). 

7T 



Using Eq. this yields 



2\/2i 



-Q(x,t) = 6+{x,t) 



2-y/2. 



-It (25) 



where 0+{x,t) describes equilibrium fluctuations (I = 
0). The induced voltage along the strip is then given by 
V = (V(L/2,t)), where the voltage operator V(x,t) is 
dictated by the Josephson relation 



V = 



V2e 



'+ 1 



(26) 



<^>+ = i[H,(j>+\. 
Using H = H Q + H D [Eqs. 0,((23j] we find 

-f E I ^'C«K) sin [20 n (a/,t)] . 



(27) 



where </> + (i) is the operator in the interaction represen- 
tation 



<p + (t) = e iHot <p+e- iHot 



and 



u(t) = e KHo+H D )t e -m t_ 



(29) 



(30) 



Assuming a weak disorder which allows a perturbative 
treatment of Hjj, u(t) is given to first order by 



u(t) = l + i J dt'H D (t'). 



(31) 



Substituting Eq. (31 1 in Eq. (28), one obtains 



+ (x,t)) =i dt'l H D (t'),(f> + (x,t) 



(32) 



Using Eqs. p3]>, |26|) and (|32|, and recalling Eq. (j6j, 
we obtain an expression for the d.c. voltage 



The time-evolution of </> + (x, t) can be expressed as 

<p + (t) = u{t)4> + (t)u\t), (28) where 



V = V! + V 2 



(33) 



^1(2 



iDLtt 







f dt '( 


sin ^ 


1 —00 





(34) 



(here 9±{t) = 0±(0,t)). Introducing the operators 

A m (x,t) ee e <VI(M<M)±Mx,t)) (35) 



where 9 + is defined in Eq. (25), we obtain the voltage- 
current characteristic 



V ^ = ^ E / ^-0{e'-< i, - i )([i n W,ii(t')])--'- (t '- t »([4(t),A n (0])} • (36) 

e 71=1,2"' °° 



r 



In terms of the retarded Green's functions we finally obtain 

*&>(t) =-ie(t)([A,(t) 4,(0)]) (37) U(J) =^ J2 Jdtsxn(?g)esm{xn(t)} 

= -2e{t)^m{Xn{t)} n=l,2 " ( 39 ) 

with =^ E J dt S m(^)Xn(t) 

n=l,2 —00 

Xnif) = ( j 4n(^)^n(0)), (38) where in the last step we have used the fact that 



G 



^sm{xn(t)} is the antisymmetric part of Xn(t) under 
t -> — t. This correlation function can be evaluated uti- 
lizing the low-energy theory developed in Sec. [TTJ 

To leading order in the perturbation Hp, the expec- 
tation value () may be replaced by ()o, evaluated with 
respect to Hq. Since the 9+, 6L degrees of freedom are 
decoupled in H , the correlation function 



Xi = X2 = X 



(40) 



where 



x (t) = / e iV2(e + (x,t)+9_(x,t)) e -iV2(e + (o,o)+0_(o,o))\ 



decouples into 



/ o 
(41) 



X(t) = Xc+(t)xc-(t) + Xs+(t)xs-(t) + Xs+(t)xc-(t)+Xc+(t)xs-(t) 

XC ± = (cos{^±(i)}cos{^±(0)}) ± , X s± = (sin{V20 ± (i)}sin{V20 ± (O)}}± 



(42) 
(43) 



Here ()± are evaluated with respect to H±. The sym- 
metric mode described by H + is a Luttinger liquid [see 
Eq. @], henceP 



Xc+(t) = XS+(t) 



lim 

e-s-0 



(iraT/v + ) 



sinh{7rT (t - ie)} 



(44) 



In contrast, as discussed below, the correlations char- 
acterizing the antisymmetric mode \xc— (t) an d XS-(t)] 
depend crucially on the parameters of (15), and in par- 
ticular on the magnitude and sign of the masses A u ^ d . 

To evaluate Xc— an d Xs- ; we nrs t note that in terms 
of the field 9 [Eq. (10)], they correspond to correlation 
functions of cos 9, sin (9, which lack a local representation 
in terms of Fermion fields. However, a convenient ex- 
pression is available in terms of the two s pecies of order 
(c-m) an d disorder (<7«,d) Ising field d 17|19 4 for A d > 0, 



cos( 



smf 



(45) 



For A d < 0, the roles of <r dl a d are simply interchanged. 
The correlators Xc--, Xs- can therefore be expressed in 
terms of C x (t) = {a x (t)a x (0)), C x (t) = (* x (t)a x (0)) 



(A = u,d), which have known analyti c approx imations 
in the semi-classical regime (|A^| S> T) 17|20 | 21 j 

C x (t)~ lA^KoiilA^t), C x (t)~ IAaI 1 / 4 (46) 

[with K (z) the modified Bessel function]. In the quan- 
tum critical regime (|A d | < T), C d (t) ~ - t" 1 / 4 . 
Employing Eqs. (44), (46), it is possible to evaluate 



the retarded correlation function and thus V{I) in either 
side of the quantum critical point of the Ising model d. 
Below we show that the resulting dramatically distinct 
behavior of the dissipative transport in the disordered 
and ordered phases of the Ising system identifies them as 
"superconducting" and "insulating" , respectively. 

A. Superconducting phases 

We first derive expressions for the I — V character- 
istics near commensurate fields Bpj [Eq. J5II where 

A d 



A { ° ] > 0, in the low T regime where Eq. <|46|) 



holds. Neglecting terms of order e 
first order in D 7 we obtain form Eq. 



'/ T and keeping the 
{391 



V {1) {I) = C 



at sin I — — 1 



(iraT/v^ 



2e / Vsinh(7rT(< - ie)) 



K (iA d t) , where C cx DL\A u A d 



1/4 



(47) 



For < Ad, this yields a non-linear I — V curve 



V4 (1) oc D[A d (B)] 



K+(B) 



(48) 



which exhibits a threshold at a critical current I c = 2eAd 
in the limit T — > 0. In the Ohmic regime I/e <C T, one 
obtains a contribution to the magnetoresistance of the 



form 



RW(T,B) w R s 
R s <xD(A d (B)) 



Ad(B) „-A d (B)/T 



(49) 



Superimposed on a moderate monotonic increase with 
B arising from K + (B) due to the suppression of p s [Eq. 
Q], the exponential factor leads to a strong decrease and 
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FIG. 2: (color online) Differential resistance vs. current bias 
in the superconducting phase for temperatures T = 0.1K, 
0.2K, 0.4K, 0.6K, 0.8K and 1.0K, for a fixed B such that 
Ad = 1.0K and K+ — 2.1 (see text); the disorder parameter 
is chosen such that Da 3 /vL = 0.1. 



flW -> at T -> as long as A^B) > is finite. The 
disordered Ising phase is thus identified as superconduct- 
ing: it corresponds to a state where the phase of the SC 
order-parameter in the two wires is locked. This suggests 
that the fields o~d physically represent phase-slips in the 
antisymmetric sector (which are gapped in this regime). 

The above analysis indicates that the first order in D 
yields an exponentially small voltage for I, T — > 0, sug- 
gesting that one should examine the perturbation scheme 
in Hjj [Eq. ( 23 )] more carefully^. Indeed, if we evalu- 



ate the expectation value ( ) expanding to the next order 
in D, we find that the correlation functions \n acquire 
corrections to Eq. ( 41 1 of the form 



S X oc fl (e±^i(t) e ±«9 ! (t) eT i2« 1 {o) eT i2« 2 (o) >o 
Using Eq.([6]), this can be written as 

The resulting contribution to the voltage 



(50) 



(51) 



DLn 

4e 



dt sin 



wit 

~2e 



5x(t) (52) 



is associated with scattering processes which do not in- 
volve the antisymmetric mode, and hence are not affected 
by the superconducting order. These correspond to co- 
incidental events incorporating two scatterers located on 
two different wires simultaneously, and therefore their 
probability is of the order of D 2 . The gapless symmetric 
mode experiences backscattering in such events, similarly 
to the usual plasmon mode in a strictly ID SC wire. In- 
serting the Luttinger liquid correlation function 



Sx(t) 

wo obtairPl 



D 



— (■KaT/v + ) 
sinh{7rT (t-ie)} 



(53) 



V^{I) = VP\B[-^- + 



il 



AeT 



2 



,1- 



- B ( — 



V 4eT 



2 4 



where K (2) oc D 2 



(54) 



r(z)rfo) 



is the Beta function. 



and B(x, y) - T{x+y) 

The full I—V characteristic in the SC phases (A l j(B) > 
T,I/e) can finally be expressed as 



V(I) = V (1) {I) + V {2 \l) , 



(55) 



where V^{T), V^(I) represent contributions from odd 
and even orders in the disorder parameter D respectively, 
and can be viewed as two resistors connected in series. 
To leading order in D, they are given by Eqs. (48) and 
(54), yielding the /, T-dependence depicted in Fig. 2. 



Note that although the second term is higher order in the 
scattering rate D, it becomes the dominant contribution 
in the limits T, / — ¥ as the first term is exponentially 
suppressed. For I/e ^> T, this indicates a power-law 
I — V relation 



V{I) - D 2 I K{B ^ +1 , K (B) 



K+{B) 



(56) 



and in the Ohmic regime (I/e <C T) 
R(T, B) - D 2 T K{B) 



(57) 



By definition of the Luttinger parameters K± [Eq. Q], 
K + > K_ and hence the assumption K_ — 2 implies 
K + > 2. As a consequence, the exponent k(B) [Eq. 
( 56 1] is small and slightly negative. We therefore con- 



clude that in spite of the phase-locking ordering of the 
antisymmetric phase mode, the true T, / — > behavior of 
the electric transport exhibits an insulating behavior. In 
practice, however, the insulating character may be mani- 
fested only at extremely low T. At moderately low T, the 
sub-leading term V^\l) is expected to be appreciable, 
and indicate a threshold at a critical current I c , directly 
related to an activation gap in the Ohmic resistance [Eq. 
(@9}] : 



log R ~ A d 



-kIc 
~2e 



(58) 
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The oscillatory nature of Ad(B) as B is tuned through 
commensurate and incommensurate values should be re- 
flected in the -B-dependence of I c , which is maximized at 
commensurate values Bn and vanishes in the vicinity of 
incommensurate regimes B ~ B N+ 1 . 



are in the ordered phase, and for T <C | A^| the correlation 
function characterizing the antisymmetric mode is given 
up to exponentially small corrections by a constant 



X-(t) ~ \ A u A d 



1/4 



(59) 



B. Insulating phases 

We next consider the insulating phase, realized in the 
vicinity of incommensurate fields B ~ B N+ i such that 
Ad < 0. In this case, both species of Ising models u and d 



As a result, x(t) = x+(t)X-(t) [Eq. (41 )] is dominated 
by the Luttinger liquid correlations [Eq. ([44 )] character- 
izing the symmetric mode. Keeping the leading order in 
D in Eq.(39), we thus find an expression for the I —V 



characteristics of the form 



V(I) = Vi [b (- 



if 
" 4eT 



1 

2 if 4 



K 4 



B 



u 

4eT 



1 

2A' 4 



A)} 



where V, oc D\A U A, 



11/4 



(60) 




FIG. 3: (color online) Differential resistance vs. current 
bias in the insulating phase for temperatures T = 0.1K, 
0.2K, 0.4K, 0.6K, 0.8K and 1.0K, for a fixed B such that 
Ad| = 1.0K and K+ — 2.1 (see text); the disorder parameter 
is chosen such that Da 3 fvL = 0.1. 



Typical plots of the resulting dynamic resistance dV/dl 
vs. / are depicted in Fig. 3, indicating a zero-bias peak 
at I —¥ 0, in sharp distinction from the SC phase (Fig. 
1). For I/e ^> T, we obtain a diverging power-law 



V(I) ~ DI 1 -^ , -y(B) = 2- 

and in the Ohmic regime (- <C T) 

R(T,B) - DT- l{B) . 



1 



K+(B) 



(61) 



(62) 



Compared to the power-law contributions to dissipation 
in the SC phase [Eqs. (56) and (57)], these results indi- 
cate a stronger divergence at low T and /. This behavior 
stems from the fact that the antisymmetric mode is in 
the insulating, charge-ordered phase, and consequently 



backscattering processes by a single impurity are favored. 
Moreover, since K+ 1 < h, the exponent j(B) > 3/2 in- 
dicating that the disorder potential is highly relevant. In 
the truly T, I — > limit (i.e., below a crossover temper- 
ature scale TJ oc which depends on the disorder strength 
D) , the perturbative treatment of Hp leading to Eq. ( 60 ) 



is not valid and localization takes over, yielding an expo- 
nentially diverging resistance^. We note that at moder- 
ately low T and /, Eq. (60) is still valid and appears to 
be compatible with the experimental data^l 



C. Critical regime 

The above analysis implies that the quantum critical 
points at Be (where = 0) correspond to SC-I and I- 
SC transitions alternately, associated with the change of 
ordering in the antisymmetric mode from phase-ordered 
to charge-ordered ground state. These transitions are 
marked by a dramatic qualitative change in the shape of 
the non-linear I — V curves, and in the T-dependence of 
the Ohmic resistance, as B crosses B^ N \ However, note 
that unlike the 2D SIT, the quantum critical points can 
not be easily identified in the transport properties, e.g. 
as crossing points of isotherms where R(B 7 T) exhibits 
a metallic behavior. In the critical regime (T 3> |A^|), 
the antisymmetric mode is characterized by power-law 
correlations %-(i) ~ t -1 / 4 and consequently 



R(T,B) - Ti-^ 



(63) 



This reflects once again an insulating behavior, charac- 
teristic to the ID nature of the system. It stems from 
the presence of a gapless mode (the symmetric plasmon) , 
which is not immune to backscattering processes. 
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IV. DISCUSSION 

In this study, we have shown that the low-T trans- 
port properties of a ladder-like superconducting device 
subject to a perpendicular magnetic field may signify a 
multitude of quantum phase transitions from a SC to 
insulating phases alternately, when its parameters are 
tuned close to the 2D SIT. These transitions stem from 
the quantum mechanical nature of the vortex chain ac- 
commodated along the central axis of the device, and 
reflect the competition between a Josephson coupling 
and a charging energy between the SC edges of the 
device, which govern the antisymmetric phase-charge 
mode. The former dominates near commensurate values 
of the vortex density, and the latter near incommensu- 
rate (1/2-integer) densities. The quantum critical points 
are of the Ising type: this is a manifestation of the Zi 
symmetry characterizing the antisymmetric mode, asso- 
ciated with interchanging the two legs of the ladder. 

The analysis presented in Sec. |III| indicates, however, 
that the electric transport properties are complicated by 
the presence of a gap less symmetric phase-charge mode, 
which provides a dissipative environment. As a result, 
the voltage response to a current bias does not exhibit 
a strictly superconducting behavior even in the phases 
classified as SC. Nevertheless, for weakly disordered sys- 
tems it is possible to observe a clear signature of the SC 
nature of these phases at finite T and /. Subtracting the 
contribution of backscattering exclusive to the symmet- 
ric mode, which can be viewed as a resistor connected 
in series, one obtains an activated behavior of the I — V 
curve and the T-dependent resistance [see Fig. 1 and 
Eq. (58)]. This behavior is sharply distinct from the in- 
sulating phases, where the differential resistance dV/dl 
exhibit a zero-bias anomaly peak [see Fig. 2]. Moreover, 
in principle it is possible to detect the quantum critical 
points (B^) separating the two phases by probing the 
B-dependence of the activated gap [Eq. (58)]. 

It should be noted that the analysis thus far relies on 
some crucial simplifying assumptions. In particular, it 
has been assumed that the model for the antisymmetric 
mode is tuned to a self-dual point, where K_ — 2. In this 
special point, where both the phase and charge fields are 
not well-defined, the chain of vortices is exactly describ- 
able in terms of free Fermions. The question arises, to 
what extent our results are robust against a finite detun- 
ing away from the self-dual point, i.e. when K- = 2+SK. 



Such corrections induce interactions among the Fermions. 
However, since in both the SC and insulating phases 
the Fermions are massive and excitations are gapped, 
these interactions can be treated perturbatively as long 
as (v—/a)SK -C |A<j|. This approximation fails when 
| Ad | — > and the critical point is shifted, but the Ising- 
type nature of the transition is maintained 25 . The phe- 
nomenology manifested by the transport properties as 
discussed above would therefore be essentially the same. 

Another point of concern when adapting the model to 
describe a realistic system is the role of finite size effects. 
In Sec. |III| the correlation functions were evaluated for 
finite T assuming that the length of the system L — > oo. 
However, we note that the SC nanowires studied, e.g., in 
Ref. 13 typically have a finite length of the order of a few 
microns. This introduces an additional low-energy cutoff 
Tl = v_/L, Using typical values of the plasma velocity 
for v- (see, e.g., Ref. [5]), we estimate Tl ~ IK. This 
implies that for sub-Kelvin temperatures, Tl effectively 
replaces T as the low-energy cutoff. In the SC phases, 
the activated contribution to the resistance is therefore 



expected to be ~ e 



-A d /T L 



Noting that Tl is also asso- 



ciated with the zero-point energy of phase-fluctuations, 
this represents contribution due to macroscopic quan- 
tum tunneling of vortices out of a metastable state in 
the finite-size SC device^. 

Finally, we wish to point out that a ladder-like SC de- 
vice where the parameters are conveniently tunable (e.g., 
a Josephson ladder) can serve as an interesting play- 
ground for the study of emergent fractional degrees of 
freedom. In particular, when the gap A^ vanishes, the 
eigenstates of Eq. ( 15 ) (at zero energy) become Majorana 
Fermions. Therefore, as recently proposed by TsvelilP^, 
inhomogeneous SC devices can be potentially utilized to 
realize localized Majorana modes at interfaces between 
superconducting and insulating segments. 
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